clc;
clear;
f=@(x)exp(-(x(1))^(1/3))*sin(sqrt(x(2)))+exp(-x(3)^2*x(4)^3*x(5));
tic,
I=quadndg(f,[0 0 0 0 0],[5,4,1,2,3]),
toc